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Abstract 

We discuss an 0{N) exension of the Sine-Gordon (S-G)equation described by the Lagrangian 



^ cos /3y^, where p — ^ which allows us to perform an expansion around the 

O ; 

. leading order in large-N result using Path- Integral methods. In leading order we show our methods 

^ ' agree with the results of a variational calculation at large-N. We discuss the striking differences 

"t^ ; 

Q^' for a non-polynomial interaction between the form for (^[0]) in the Gaussian approximation that 

. one obtains at large-N when compared to the = 1 case. This is in contrast to the case when 

, is a polynomial and no such drastic differences occur. We find for our large-N extension of 

' the Sine-Gordon model that the unbroken ground state is unstable as one increases the coupling 

! 

constant (as it is for the original S-G equation) and we find in leading order that the unbroken 
symmetry vacuum is stable as long as < 247r. 
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I. INTRODUCTION 



Although the Sine-Gordon equation has previously been studied using variational meth- 



ods, 



a systematic expansion around the mean field theory result has been lacking. 
Here we use the method of auxiliary fields to study going beyond mean field theory using a 
systematic expansion based on the the parameter obtained when we consider a particular 
0(N) extension of the orginal S-G model. Using the methodology of P| we first extend 
the problem to having 0(N) symmetry and then introduce an auxiliary field p = by 
the insertion of a functional delta function into the path integral. The large-N expansion is 
then obtained by integrating out exactly the (p held and then performing a steepest descent 
evaluation of the remaining path integral. Unlike the case of polynomial interactions, at 
leading order in the Gaussian approximation there is a distinct difference in the equations 
of motion when = 1 and A^ large. We display this by comparing the 0^ results with those 
of the S-G equation. As in the original S-G equation, our 0{N) extended model has an 
unstable vacuum as a function of the coupling constant. We find that in leading order in 
large-N that the vacuum is unstable for > 247r. 

II. TIME DEPENDENT VARIATIONAL APPROACH 

One method for obtaining a time dependent variational approximation to a quantum 
field theory is to start with Dirac's Variational principle pfor obtaining the functional 
Schrodinger Equation. By using wave functionals which become "exact" in the large-N 
limit- namely Gaussians we will obtain the time- dependent Hartree approximation to the 
exact field equations. 

Dirac's variational principle 

r = J dt< -i^i^ > (1) 

leads directly to the Schrodinger equation: 

ST = 0^{t^^-H}\^>=0 (2) 

In the (fi representation 



2 



and 



H^J rx[~S'/[Scp{x)r + ^Vcp{x)Vcp{x) + VM (3) 
In this representation one can choose a Gaussian trial wave functionah 

< ifil'^y >= il'v[(p,t] = Aexp{- / {[(p{x) - 4>c{x,t)] 

Jx,y 

[G-\x,y,t)/A-i^x,y,t)My)-My,t)]} + i 1 7rc(x, - i)]}, (4) 

Jx 

where (f)c{x,t) — {^v\(p{x)\^y); ndx^t) — (^^| — iS/S(p{x)\^y) and 

G{x, y, t) = {^^\ip{x)ip{y)\^^) - (P^x, t)(j)c{y, t). (5) 
Then the effective action for the trial wave functional is 

r(0c,7rc,G',E) = j dt<'^^\id/dt-H\^!^> 

— j dtdx['Kc{x,t)d(f)c{x,t)/dt + j dtdxdy'E{x,y;t)dG{x,y;t)/dt] 

- J dt<H> (6) 

where 

<H> = f dxdir^ + 2[J:Gmx,x) + lG-\x,x) + l-{V(p)^ 
J 2, 8 2 

+ ^limV.V,G(a;,2/)+< V>}. (7) 

< if > is a constant of the motion and is a first integral of the motion. The equations 
one gets by varying the effective action with respect to the variational parameters are: 

2 , . d<V> 



0(Pc 

4)c{x,t) =7:c{x,t) 
G(x, y;t) = 2 1 dz[J:(x, z)G(z, y) + G{x, z)T,{z, y)\ 

t{x,y;t) = J dz[-2J:{x,z)E{z,y) + ^G-\x,z)G-\z,y)] 
For calculating expectation values, 

p[0] = >i^exp{- / [^{x) - M^M^^^^^^iviv) - MvM (8) 



with the normahzation A determined from the conservation of probabihty: 

J Vct> P[<P] = 1 (9) 
To get a general expression for (V) one can use the method of the generating functionah 



Z[j] = {expll j{x,t){^{x,t)-4>,{x,t)}])=exp[^ / j{x,t)G{xy,t)j{y,t)]. (10) 



In the single field case (N=l) for example one has 

6j'^^{x 



((0(x,t)-0e(x,t))^") = -^|,=O. (11) 



and we obtain the simple result 

mx, t) - Ux, t)f^) = (2n - l)!!G(a;, x- (12) 

To obtain the expectation value of an arbitrary potential we assume that the potential has a 
Taylor series expansion about the classical field 0c = (0) and the expectation value is taken 
with respect to the Gaussian trial wave function. At N=l we obtain the result: 

^^^"^^^ - \ {2n)\ dct^l- ^^"^-"^^^ ^^^^ 

For the 0(N) extension of the potential, (p ^ (pi] i = 1 . . . N and the Green's function G is 
now an X matrix. One then finds making the Gaussian ansatz and ignoring non-leading 
terms in 1/A^ that (see j?!) 



where 



{{^r) - ip)" (14) 



W-^V^-G+^-G + e (15) 



Thus at leading order in large- A^ the expectation value of the potential becomes 

{V[^]) = VM + G] (16) 

A. (j)^ Model 

First let us look at a typical polynomial model such as the 0^ model described by the 
classical potential 



For N — 1 one gets the Hartree approximation result 

{V) =^[02 + G] + ^[0^ + 602G + 3G^] 

+?[0^ + 150^G' + 450'G'' + 15G''] (18) 
6 

where in the above G — G{x, x) and G{x, y; t) is the variational function. 

To obtain the phase structure one calculates the effective potential which is the expec- 
tation value of the Hamiltonian in the Gaussian trial state for spatially homoogenous time 
independent fields: 

Veff = + ^lim,^yV,VyG{x,y) + {V). (19) 

We can change variables from G{x, y) to the effective mass x for constant (homogeneous) 
external sources and in that case parametrize G{x,y) via 

G(..,)./ge<'<-)5^. (20) 



So that if we introduce a cutoff A 



1 

G(x,x)^—ln— (21) 
47r X 



The first two terms in the effective potential coming from the kinetic energy can now be 
written as 

This has an infinite x independent term which is the cosmological constant and needs to be 
subtracted by hand. 

In general Vgff is a functional of 0, G or 0, x- Using the chain rule we have 

dV ^ dVdG 

dx dGdx ^ ^ 



with 
Here 

dV 1 



- - (24) 
dx 47rx ^ ^ 



^[-X + ^l' + m<l>' + G[x]} 

+5r7{(/.^ + 6</.^G[x] + 3G2[x]}] (25) 



Another way of writing the Kinetic energy which removes the cosmological constant term 
is used in Q and leads to a result independent of the cutoff. There one lets x ~^ i^i fhe 
definition of G and sets 

The potential energy also has divergences which are related to mass and coupling constant 
renormalization if one does not do any normal ordering. Here we will first use conventional 
renormalization and then show it is equivalent (in 1+1 dimensions) to normal ordering the 
unrenormalized result with respect to the renormalized mass. First let us look at the case 

= 1. We define the renormalized mass as the value of x determined from the gap equation 



obtained from setting l^|,/,=o = i.e. 



m 



2 



ix^ + SAoGim^] + lhriG^[m^] (27) 



The gap equation for x is fully renormalized if we solve the above for n in terms of m and 
also perform the coupling constant renormalization 

= Ao + lQr]G[m^] (28) 

Defining 

GJ^] = GW]-G\m^] = 

4vr X 



1 IT) ^ 

GM = G[x]-G[m'] = -\n— (29) 



one then obtains the renormalized equation 

1^ =\[-X + ^' + ^K{ct>' + Gr[x]] 

+5ry{0^ + 60UM + 3G^[x]}] (30) 

We notice here by comparing Eq. ESI and Eq. EHl that the only difference between the 
unrenormalized and renormalized gap equations are the replacement of: 

fi —> m; X^Xr;G^Gr- (31) 

n 

This is an example of the result in two space-time dimensions (see also 8]) that if one did 
no mass and coupling constant renormalization, but normal ordered the two point function 
with respect to the renormalized mass parameter m the theory would be rendered finite. 

To obtain the renormalized effective potential one needs to do an integration of the 
renormalized equaition for ^ and add the pure dependent terms. 



We have 



dV 
dx 



■[-X + m^ + 3Xr{<P' + Gr\x\} 



+5n{cf>^ + 6cf>'Gr[x] + SGl[x]}] 



Integrating and adding the pure 4> dependent terms we obtain: 

z 4 o 

1 771 ^ 

— [x-m^ + (m^ + 3A,02 + brjcj)^) In — ] 
Stt X 

^(A. + 10#^)ln^- + ^,ln3_. 



(32) 



(33) 



B. ^^at large AT 

At leading order in large N one has as discussed before that 

{V)/N = ^f,'[cl>l + G{x, x)] + ^[cl>l + G{x, x)r + + G{x, x)]\ (34) 



where here we use the shorthand 0^ 



now: 



AT 



. Performing the mass renormalization we find 



— ji^ -\- XoG{x, x; m^) + r]oG'^{x, x; w?). (35) 
The equation for the coupling constant renormalization is also slightly changed and reads: 

Ar = Ao + 2r]G{x, x; w?) (36) 



[-X + + \r{<f + Gr{x, x] x)] + ^{4>^ + Gr{x, X] x)Y 



(37) 



then the renormalized gap equation is 

9y _ 1 

dG~T 

Thus apart from the factors in front of A and 77 reflecting 3 = 1+2 /N etc, the renormalization 
at large N is quite similar to what happens in the Hartree approximation. We also find 

dH dH dG 



dx dG dx 

which after renormalization can be written (using for i7 at constant fields) 



(38) 



eff 



dx Sttx 



X-m^ - Xr 



(.2 + ]-r][<p' + -^ln[— 



2 \ 2' 



(39) 



This result again could have been obtained if we did NOT renormalize the coupling con- 
stant or the mass and instead used m? for the mass parameter and everywhere made the 
subtraction: 

G{xx^ X) G{xx, x) ^ G{xx, rn?) = Gr{x, x; x)- (40) 
Integrating as before we find the renormalized effective potential can be written in the form 

V[<P,x]/N= ^[x-m'] + lm'[<p' + Gr[x]] 

+iA,[02 + G,]2 + W + af (41) 

4 D 

which has the advertised form. 



C. sine-Gordon equation 

The sine-Gordon equation is quite interesting since the classical equation has a kink 



solution which remains after quantization and there is a. 
the effective coupling constant as discussed by Coleman 
is described by the Hamiltonian: 



so a phase transition as we increase 
. The usual sine-Gordon equation 



H = |rfx[^7r(x,t)2 + i(Vy.(x,t))2-pcos/5y. + 7] (42) 
where in the "93" representation 

Taking the expectation value of this Hamiltonian in the trial Gaussian wave functional 
state one now has 

< V >=< cos/?v7 >= cos/50 e-^^("'") (43) 

Varying the action we obtain the following field equations for the variational functions: 

0- V20+ ^sin/30 e-^^("'") = 0, (44) 

G(x, y) = 2j dz{G{x, z)^{z, y) + S(x, z)G{z, y)} (45) 

t{x,y-t) = J dz[-2J:{x,z)i:{z,y) + ^G-\x,z)G-\z,y)] 

+ [^V^-^«ocos/50(x,t) e-^^(^'^)]5(x-y) (46) 



In the Heisenberg picture, making the same the Hartree approximation gives the following 
covariant equation for the nonequaltime correlation function: 



P + {^)]g{x,x')=6\x-x') (47) 



where 

K2 



(^) = «ocos/30(x,t) e-^^(^'^) = x{x). (48) 



Again the connection between the covariant and noncovariant Green's functions at equal 
spatio-temporal points is: = G{x,x;t). 

When X is independent of x one can again introduce the Fourier Transform: 

G{x,y) = ^ I dkG{k)e'^^''-y^ (49) 

ZTT J 

where now 

G{k) ^ 



2(A:2 + x)i/2 

In particular, when = one has the gap equation: 

^ = ao e-'^^^^'^'^^) = ao(— )-^ (50) 

This equation tells one how to choose as a function of the cutoff A to insure that the 
physical mass x is finite. 

To study the phase transition in this theory in this approximation we follow j^. We 
determine for what values of the coupling constant one can have = be a minimum of 
the energy. 

In the vacuum sector we have for the Hamiltonian Density: 



<H>= \g-\xx- x) + i lim V.V,G(x, x) - ^ e'^^^^'^^^) (51) 
Inserting a momentum space cutoff A 



i-nJ-K (A;2+ ^2)1/2 ^2V^2 

Let us define an arbitrary renormalized mass squared parameter by 

ar{m') = aoi^f^'- (53) 

which has the property that 

«r(x) = X (54) 



we find that the once subtracted (at m^) energy density is given by 

nix) - n{m^) = ^ix - m^) - %{^ri'- (55) 



The first derivative is given by: 

dx Stt 
The minimum for </> = is at 



(56) 



(57) 



We notice that if we choose 

= ar(m^) 

then the extrema of the energy is at 

X = m?. 

Since the second derivative of the energy is given by: 

d^X 8nm^^87i '^m^' ^ ' 

we find that a stable ground state with unbroken symmetry (0 = 0) exists only for 

!— < 1 (59) 

which is in accord with the result of Coleman Q|. 
The resulting renormalized effective potential is 

V[<P,X] = ^Ix-m'] - ^{^Y^'/'-^cosP<l> (60) 

For the Schrodinger picture update equations we can defined a finite spatially and tem- 
porally varying renormalized mass M^(x, t) via 

M^(x, t) = exp[— — {(^(xx; x) — G{xx; m^)}] = exp[— — (^^(xx; x)] (61) 

This mass renormalization will render the resulting update equations for the time evo- 
lution problem finite. In terms of the renormalized mass M^(x, t) the renormalized TDHF 
equations are : 

0_V20 + ^^^!Msin/30 = o, (62) 

P 
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y) = 2j dz{G{x, z)J:{z, y) + z)G{z, y)} (63) 

t{x,y-t) = J dz[-2J:{x,z)i:{z,y) + ^G-\x,z)G-\z,y)] 

+ - ^M\x,t) cosP<P{x,t)]5'ix - y) (64) 

In the Heisenberg picture, the covariant time dependent equations in terms of M^(x,t) 

are 

[□ + M2(x,t)cos/?0(x)]^(x,x') = 5'^{x-x') (65) 
thus the space and time dependent effective mass is 

x(x) = cos (3(f){x)e~~ 2*' = M^(x, t) cos (3(f){x) (66) 

The expectation value of the field obeys the equation 

□0(x) + M^{x, t) sin (3<p{x) = 0. (67) 

The subscript r in Eq. IHSl means that one uses the once subtracted covariant Green's 
function. Here we see that the main effect of the Hartree approximation on the kink evolution 
equation is to replace by a self consistently determined M'^{x,t). We notice that the 
"classical" field evolves quite differently than the propogator. As we shall see next, making 
the Gaussian approximation but keeping only leading order terms at large-N dramatically 
alters the evolution equation in a manner that can be interpreted as a change in the behavior 
of the periodicity itself, but with the evolution of and Q being similar. For the N-component 
field we will consider instead the classical potential: 

= -N^ cos/3v/p; P^^- (68) 

Using the previous result at large that {H) = K[x] + Vd [4>l + G[x]], we have that the 
effective potential is now: 



Veff/N = ^ - ^ cosP^^l + G[x]. (69) 



The renormalized effective potental is thus given by 



Veff/N = ^-j,cos P^^l + Gr[x;m^]. (70) 
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with an arbitrary subtraction point. If we choose the extremum of energy when 
to be at X = than the renormahzed effective potential can be written as: 



^[0, X] = ^[X - -V ^ cos (/5y^0^ + I H^]^ (71) 
We notice now that the effective mass (pole in the propagator) is at 



^ sin(/3^0^ + ^ln[ig]) 
Y = m , — (72) 



4n I- X ' 

which is functionally totally different than the = 1 Hartree result that 

g'^Q^, (x X 'x) 

x(x) = cos/30(x)e 2» = M^(x, t) cos/90(x) (73) 



III. PATH INTEGRAL APPROACH 



To obtain the large-N expansion for an arbitrary even polynomial interaction we follow 
the work of Eyal et. al. j^]. Starting from the usual path integral for the generating 
functional 



Z[j] = J D<f) Exp 
one introduces the composite field p using : 

dx6{p 



NV\ 



N 



-]+j-4>} 



(74) 



dpdx exp 



(75) 



This allows us to rewrite the generating functional (also adding sources for x and p) as 
Z = J D4>DxDpExp [z{^[dM' - \x4>-$- NV[p\ + ^px + ]-$+ NSx + NJp} 

The equations resulting from the second form of the action are 

[□ + X]<P^ =J^■, X = -2J - 2V'[p] ; p = ■ 0/iV - 25 



(76) 



(77) 

The large-N expansion is obtained by integrating out the (f) fields and, recognizing that 
the resulting action is proportional to performing the remaining integral by Stationary 
Phase. The first integration yields 



Z[j,S, J]= J DxDpe'^^^ff^f''^'^'^"'^ 



(78) 
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where 

Seff = ^{[ / ^^^TT^^t^' ^] - + + Sx + Jp+ lTrLng-'[x, y, x]} 



(79) 



and here 

g-'[x,y,x]^P + xW{x-y) (80) 
is now the covariant two time Green function. The stationary phase conditions are: 

where (pc — I Gj and 4>l = 4> ■ 4>/N. In the absence of sources, this leads to the constraint 
equation for p 

p^(t>l + -Tr^? (82) 
and the equation for the effective mass (gap equation) in leading order: 

X^-2V'[p]. (83) 

At leading order we also have from the definition of 0c = y ^^^^ 

[n + x]0c = J (84) 

so that unlike the Hartree case the mass for the classical field (j) is the same as the mass 
entering into the propagator. The Gaussian fluctuations about this mean field (Stationary 
point of the effective action) is obtained from the second derivative of the effective action at 
the stationary point. This matrix is actually the inverse propagator of the two component 

field made up of x aiid p. The stability of the largc-N approximation is related to the 
eigenvalues of this matrix being positive. Critical behavior is related to the deteminant of 
the matrix being zero. The second derivatives are given by 

'i^^Mm ^ ^ "^'^^^^^^^ '^^'^'^ " ^ ^^^^ 

1 ^''S 1 5''S 1,, 

-5{x - y) (86) 



i 5x{x)6p{y) i 6x{x)6p{y) 2 
1 S^S 



i 6p{x)6p{y) 
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-V"[p]S{x-y). (87) 



For constant fields, the Fourier transform of the inverse propagator is given by 

/ S(p) 1 



2 2 



(88) 



V 2 

Expanding S about the stationary phase point, keeping up to Gaussian fluctuations and 
Legendre transforming one has the effective action to is 



F 



dxCd[(l),p,x] + i^Tr\Yig-^[x,y,x]} + -;^Tr\nD~^[x,y,x^ pA]} (89) 



and the inverse matrix propagator D~^is defined in terms of the second derivatives of the 
effective action at the stationary phase point with respect to the fields p and x- 

In 1+1 dimensions, S(g) is finite and (for the case of zero classical field) given by the 
integral 

^ = \ ^ f fPE_ 1 . V 

^> I J (2vr)2 (-p2 + _ q)2 + ^) j (2,^)2 (p| + x){pe - qE? + x) ^ ^ 

where in the last equality we have made the Wick rotation po ~^ ipoE- The value of S(0) is 
important in studying the critical points of the theory where the determinant of the quadratic 
fluctuations around the mean field becomes zero and one has a massless excitation. At that 
point the large-N expansion breaks down and a resummation technique is needed as discussed 
in ^. Explicitly 

S(g = 0,x) = ^. (91) 
and the condition for the determinant to vanish becomes 

mV"[p] + 1 = 0. (92) 



A. (p^ model at large-N 



0: 



We start with the usual Lagrangian for the 0^ model is (see j7[) 

^ = li^<l>^)' - ^/^v- - ^(0 ■ $r - ^(0 ■ 0)' (93) 

and i = 1 . . . N. 

Introducing the functional 6 function that p = ^ we get the second form of the La- 
grangian 

^ 1 /Q - ^2 1 2 ^ 2 V 3 , X, 0■0^ /o.x 
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The effective potential obtained at leading order in large N for 0^ theory written in terms 
of the auxiliary fields x, p is simply (scaling out the N) 

'elf = lpi[.(l>^ -P) + 2/^> + -fP^ + p '' ' ^ 

where here Q is the covariant two time Green function 



X(0^ -P) + ^P + 1-p' + Ip' + ^trln^-^[x] 



(95) 



g-\x,y;x) = p + x{xW{x-y). 



The gap equation is obtained from 



and p can be eliminated in favor of x using 



eff 



^ P^ (p'^ + -G{x,x;x) 



(96) 



(97) 



(98) 



dx I 

However, ^Q{x,x;x) = G{x,x;x) = ({0(^) ~ (0(^))}^) where G{x,y) was defined earlier in 
the variational approach. This is seen by considering 



1 



tS{x,x;x) 



1 r (fkE _ 1 /■ 
(kl + y) = Mi 



dk 



I • - (27r)2y {kl + x) (27r)y 2^^+^ 

Once we rewrite solely in terms of x as [0^, x, p[x]], then we obtain that 

dVeff _ d\^dp_ 
dx ^ dp dx 

This yields the unrenormalized equation: 



(99) 



(100) 



(101) 



using the same mass and coupling constant renormalization as in the variational aproach we 
get exactly the same result: 



eff 



dx Sttx 



X — m? — \r\4>^ -\- -;— ln[ 



m 



+ — ln[ 



m 



47r X / V 47r X 
Integrating this with respect to x we obtain the renormalized effective potential: 

K//[0,X] 



(102) 



-m 
2 



6 



^ [X - + (m' + A,0' + r](t)^) In ^] 
X 3847r3 X 



647r2 



(A^ + 27702) In^ 



(103) 
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We can rewrite this as: 

VefficP, X] = ^[X - m'] + \m\<\>^ + a[x]] 

+ lA,[02 + aF + lr?[02 + af (104) 

4 D 

This is exactly what we got from the variational method. However now we can systematically 
improve on the variational result in terms of a series in XjN . 

To study the critical behavior of 0^ and the stability of the large-N expansion one needs 
to study the condition for the determinant of the Gaussian fluctuations to vanish, namely: 

my"\p\ + ^ = 0. (105) 

where 

S(g = 0,x) = T^. (106) 
47rx 

Here the potential in terms of p is given by 

m-^p^^p'^\p' (107) 

The second derivative is to be evaluated at the stationary point where x — P — 

0' + ^[x, 0] 

We can rewrite everything in terms of renormalized parameters by having /i ^ m; Aq — > 
Xr and p— > (j)^ + Gr where everything is normal ordered with respect to the renormalized 
mass of the field m. Taking the second derivative at the stationary phase point we obtain 
when 0^0 and x — that (since Gr[x = m] = 0) 

V"\p] = \r + r]p^\r (108) 

. Thus the critical condition for the existence of a zero in the inverse propagator is 

+- = 0. (109) 



In (f)^ we can have the renormalized Aj. negative and still have a positive definite theory as 
long as ?7 > 0. 
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B. 0(N) Sine-Gordon 



Here the Lagrangian divided by N written in terms of the auxihary fields x, p is given by 

^[dM"-^X$-$-V{p] (110) 



where 

Following what we did for 0^ field theory we have that the effective potential is now: 



V{p] = -^cos/3Vp (111) 



Veff = lx{<l>' -P)- pcos/3v/p + ^Jrhng-'[x] (112) 
where again Q is the covariant Green's function with 

Q-' = P + x]5\x-y). (113) 
The gap equation is obtained from 

and p again can be eliminated in favor of x using 

^ = O^p = 02 + ^6;(x,x;x) (115) 

As in the polynomial potential cases, we can render the theory finite by just regulating p 
with respect to an arbitrary mass parameter mi 

PR = 4>^ + X] x) - 0{x, x; ml)] (116) 



so that 

PR 



2 + -^ln[^] (117) 



47r X 

The gap equation for % at = gives the renormalized mass m^: 

n? = ao \ , ^ - (118) 



/3^/^Mi] 

so that choosing our subtraction point to be mi = m we obtain that 



ao = (119) 
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Next we have again 

dV _dV dp 
dx dp dx 

Using our renormalization scheme p ^ pn, we obtain 



(120) 



dx 



sm 



X-m 



;i2i) 



Integrating this with respect to x gives the renormahzed effective potential 

.2 / I 1 



1 



m 



,171^ 



X] = -^[X - m'] - — cos ( + _ ln[— ] 



(122) 



/52 --^'-y- ■ 4^--- ^ 

which again displays the usual result at large-N that the renormalized expectation of the 
Classical potential has the property: 



(123) 



The second derivative of Vg// with respect to x at the minimum x = ""^^ (for the unbroken 
vacuum = 0) is 

/?2 -2471 , , 

^ (124) 



Thus the vacuum is stable as long as 



1927r2x 



P'^ < 247r. 



(125) 



Let us now show that this result is the same as demanding that the large-N expansion be 
an expansion about a minimum by considering the determinant of the fluctuations around 
the minimum. Critical behavior is given by the vanishing of the determinant (Eq. 1221) • We 
have, using our renormalization scheme, 



where now 



and pr 



m 



Vr[pr\ = -^COS/3v^ 



(126) 



(127) 



+ Gr with the subtraction being made at the physical mass m. Taking two 



derivatives we find that 



0; 



cosx 



smx 



;i28) 
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where x = (3y/p^ ^ Pr = ^lii^ At the stationary phase pont we need the hmit of V" as 
X ^ This is 

V"[p] - -m'^ (129) 
The condition for the determinant of the fluctuations to be positive is 

mV"[p] + l>0, (130) 

where 

^q = 0,x) = ^- (131) 



Eq. IIHOI leads to the condition (when x 



m 



2^ 



< 247r. (132) 

This is the same result as that we found for the stability of the vacuum 

To go beyond the lowest order approximation one first includes the Gaussian fluctuations 
and considers p and x ^is components of a new two component field with matrix inverse 
propagator D^^. The one particle irreducible generating functional now gets a contribution 

^TrLn[D-'] (133) 

However this naive 1/A^ expansion is secular for time evolution problems as shown in j^. A 
better approximation which avoids secularity is a resummed next to leading order in 1/N 
expansion obtained from the two-particle irreducible formalism [l^ written in terms of the 
fields 0, X, P and propagators Q and D. The self-consistent 1/A^ correction is found from the 
Graph in T2[^,G] 

^2 = ^Tr J dxdyg{x,y)D{x,y)g{y,x) (134) 

which will give the same Schwinger Dyson equations as the standard 1-Pl 1/N expansion 
except with FULL propagators and not leading order ones in the self energy graphs. This 
is discussed in detail in lllll m 



IV. CONCLUSIONS 



In this paper we have shown how to generalize the Sine- Gordon equation to the case of 
0(N) symmetry. This then allows us to go beyond the Gaussian approximation (at large-N) 

19 



in a systematic, controlable way. We found that at leading order in large N, the pole in 
the propagator looks functionally different from that found in the Hartree approximation for 
= 1. Thus the combinatorics of large-N for non-polynomial potentials leads to expressions 
for the evolution of the one and two point function which look qualitatively different from 
what is found for polynomial potentials. We also found that the stability of the mean-field 
vacuum found using a Hamiltonian approach leads to the same condition as the stability 
of the large-N approximation found by studying Gaussian fluctuations about the leading 
order. We briefly discussed the fact that the naive expansion is secular beyond leading 
order and has to be replaced by a resummed 1/N expansion obtainable from a two particle 
irreducibility approach. 
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